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Abstract: Recent experimental data is indicative of the existence of sterile neutrinos. 
The minimal scheme that can account for the data and is consistent with cosmological 
observations is the 3+1 picture which consists of three predominantly active and one pre- 
dominantly sterile neutrino with the fourth neutrino being heavier than the other three. 
Within this scheme there are two possibilities depending on whether the three light states 
obey normal or inverted hierarchy. In this paper we consider the two zero textures of the 
low energy neutrino mass matrix in presence of one additional sterile neutrino. We find 
that among 45 possible two zero textures for this case, 15 are consistent with all current 
observations. Remarkably, these correspond to the two-zero textures of a three active neu- 
trino mass matrix. We discuss the mass spectrum and the parameter correlations that 
we find in the various textures. We also present the effective mass governing neutrinoless 
double beta decay function of the lowest mass. 
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1. Introduction 

Neutrino oscillation in standard three flavour picture is now well established from solar, 
atmospheric, reactor and accelerator neutrinos. The mass squared differences governing 
these oscillations are ~ 10~ 4 eV 2 and 10~ 3 eV 2 . However, the reported observations of Vn — 
v e oscillations in LSND experiment Q and recent confirmation of this by the MiniBooNE 
experiment |2[ with oscillation frequency governed by a mass squared difference ~ eV 2 
cannot be accounted for in the above framework. These results motivate the introduction 
of atleast one extra neutrino of mass ~ eV to account for the three independent mass scales 
governing solar, atmospheric and LSND oscillations. LEP data on measurement of Z-line 
shape dictates that there can be only three neutrinos with standard weak interactions 
and so the fourth light neutrino, if it exists must be a Standard Model singlet or sterile. 
Recently this hypothesis garnered additional support from (i) disappearance of electron 
antineutrinos in reactor experiments with recalculated fluxes Q and (ii) deficit of electron 
neutrinos measured in the solar neutrino detectors GALLEX and SAGE using radioactive 
sources Q. The recent ICARUS results || however, did not find any evidence for the LSND 
oscillations. But this does not completely rule out the LSND parameter space and small 
active-sterile mixing still remains allowed ||. Thus, the situation with sterile neutrinos 
remains quite intriguing and many future experiments are proposed/planned to test these 
results and reach a definitive conclusion |7|. 

Addition of one extra sterile neutrino to the standard three generation picture gives 
rise to two possible mass patterns - the 2+2 and 3+1 scenarios ||. Of these, the 2+2 
schemes are strongly disfavored by the solar and atmospheric neutrino oscillation data 
|| . The 3+1 picture also suffers from some tension between observation of oscillations 
in antineutrino channel by LSND and MiniBooNE and non-observation of oscillations in 
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the neutrino channels as well as in disappearance measurements. However, it was shown 
recently in [1C] that a reasonable goodness-of-fit can still be obtained. Although introduc- 
tion of more than one sterile neutrinos may provide a better fit to the neutrino oscillation 
data [11], the 3+1 scheme is considered to be minimal and to be more consistent with the 
cosmological data [^] . Very recently combined analysis of cosmological and short baseline 



(SBL) data in the context of additional sterile neutrinos have been performed in [13, 14 



The analysis in [13] found a preference of the 3+1 scenario over 3+2 while the analysis in 
1 14] shows that the status of the 3+2 scenario depends on the cosmological data set used 
and the fitting procedure and no conclusive statement can be made regarding whether it 
is favoured or disallowed. In fact the current cosmological observations of an weakly inter- 



acting relativistic "dark radiation" may actually prefer an additional sterile neutrino [15]. 
If this radiation is attributed to extra neutrino species then the data gives the bound on 
the number of neutrinos as N e ff = 4.08 ± 0.8 at 95% C.L. |l5| . 

In this paper we consider the structure of the low energy neutrino mass matrices in 
presence of one extra sterile neutrino. The low energy mass matrix in the flavour basis is 
now complex symmetric with 10 independent entries and can be expressed in general as, 
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For three flavours the last row and the column would be absent and the mass matrix 
would contain 6 elements. In the context of three generations, a very remarkable result was 



obtained in [16] that there can be at the most two zeros in the low energy neutrino mass 
matrix in the flavour basis. For the three neutrino mass matrix there can be 15 possible 
two zero texture structures. These are the same as those shown in Table [l| for the four 
neutrino case after omitting the fourth row and column. Among these, only 7 textures 
corresponding to the A, B and C class were found to be compatible with the experimental 
data on neutrino oscillation fllql . Normal Hierarchy (NH) was found to be allowed in all 
the textures whereas Inverted Hierarchy (IH) and Quasi-Degenerate (QD) solutions were 
allowed for the B and C classes. Various aspects of the two zero textures in the low energy 



neutrino mass matrix have been examined in [17, 18]. Recently, this analysis has been 
redone in IS, 21, 22, |23|] to take into account the recent results including the measured 
values of mixing angle #13 by the reactor experiments [p4f |. The analysis including the 
latest data and allowing the parameters to vary randomly in their 3<r range shows that all 
7 textures of the original analysis in [16| remain allowed p3[j . However, with the la range 
of parameters the scenarios become more constrained. With the oscillation parameters 



taken from [25] only A class with NH remain allowed while for oscillation parameters in 



[26] the textures belonging to B2 and -B4 classes for IH and C for NH get excluded [23]. 
This demonstrates that with precise determination of oscillation parameters the allowed 
scenarios would become more constrained. However, the situation may change altogether 
in presence of sterile neutrinos. 
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In this paper, we examine how many two zero textures are allowed by the current 
oscillation data in the low energy neutrino mass matrix when an extra sterile neutrino 
is present. We assume the known oscillation parameters are normally distributed with 
the peak at the best-fit value and la error as the width. First we check the status of 
the two zero texture solutions in the context of three generation mass matrices by this 
procedure. Then we check how much these conclusions change in the 3+1 scenario with 
one additional sterile neutrino. We also investigate if any new interesting correlations can 
be found specially for the sterile mixing angles. Finally, we discuss the changes expected in 
our result if the mass and mixing parameters are varied randomly in their 3a range instead 
of varying them in a Gaussian distribution peaked at the best-fit value. 

Texture zero implies some of the elements are much smaller than the other elements 
of the mass matrix. Analysis of texture zeros puts restriction on the nature of the mass 
spectrum and can give rise to correlations between the mixing angles, masses and CP phases 
which may be confirmed or falsified by experimental observations. This can often help in 
understanding the underlying flavour symmetry |27j. In case neutrino mass is generated 
by seesaw mechanism the texture zeros in the low energy mass matrix can be useful in 
identifying the possible high scale Yukawa matrices [28, 29, |3C|], 



The plan of the paper goes as follows. In the next section we discuss our formalism. 
Section III discusses the results. We end in section IV with summary and conclusions. 



2. Formalism 

We consider the 3+1 mass spectrum. This can generate two possible mass orderings. In 
one case the fourth neutrino is heavier than the other three and in the other case the 
fourth neutrino is lighter than the other three. LSND/MiniBooNE observations dictate 
that the mass squared difference of the fourth state with the three other states is ~ eV 2 . 
However, the scheme in which the fourth state is lower would be more disfavored from 
cosmology since there will be three neutrino states with mass ~ eV which will contribute 
to the cosmological energy density. Therefore, we consider the picture in which the fourth 
state is heavier. Then, there are two possibilities shown in Fig. 1. 



1. mi ~ m 2 < m 3 < m A which corresponds to a normal ordering among the active 
neutrinos (SNH). This gives 

m 2 = ^/m 2 + Am 2 2 , m 3 = ^m 2 + Am 2 12 + Am 2 23 , m 4 = ^m 2 + Am 2 4 . 

2. 7773 < m2 ~ m \ < m 4 corresponding an inverted ordering among the active neutrinos 
(SIH) with the masses 

771-1 = J 17^3 + ^ m 13 > m 2 = \j m \ + A m 13 + Am 2 2 , TTI4 = J 771 2 + Am 2 ^. 

Here, Am?- = m 2 — mf. 

We assume that the charged lepton mass matrix is diagonal and the mixing in the neu- 
trino sector is solely responsible for the leptonic mixing. In the present case, the neutrino 
mixing matrix, V can be parametrized in terms of six mixing angles (#13, #12, #14,#23,#24)#34), 
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Table 1: Allowed two zero textures in the 3+1 scenario. The 15 possible two zero textures of three 
neutrino case are same as these after omitting the 4th row and column. 



three Dirac type CP phases (813, $14, #24) and three Majorana type CP phases (a,/3,7). The 
neutrino mass matrix in flavour basis is given by 

M v = V*Mf a9 V ] (2.1) 

where, M^ mg = Diag{m\, rri2, 1713, 7714). 



V = U.P H with 

U = -R34-R24-R14-R23-R13-R12 (2.2) 
where Rij represent rotation in the ij generation space, for instance: 



- 4 - 



Am 



Am 



Am 



SNH 



Am" 



Am" 



Am" 



SIH 



Figure 1: The allowed 3+1 mass ordering 
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with Sij = sinOij and Cy = cosOij. The diagonal phase matrix has the form 



P = Diag{\, e - fa /2 > e -i(^/2-*i 3 ) j e -<( 7 /2-«w)). 

The best-fit values and the lcr and 3cr ranges of the oscillation parameters in the 3+1 
scenario are given in Table 2. One can define three mass ratios 



m,2 rris mi 

The two zero textures in the neutrino mass matrix give two complex equations viz. 

M u{ab) = 0, 
M u(pq) = 0. 



(2.3) 



(2.4) 



where a, b, p and q can take the values e, r and s. The above eqn.(2.4) can be written 
as 

U al U bl + -U a2 U b2 + -U a3 U b3 e 2iS ^ + zU a4 U M = 0, (2.5) 
x y 
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Table 2: The experimental constraints on neutrino parameters. The three generation constraints 
are from global analysis in p5[ . The constraints on sterile parameters involving the fourth neutrino 
are from . The constraint on sin 2 6*34 is obtained from |32j . 



U pl U ql + -U p2 U q2 + -U p3 U q3 e 2iS ™ + zU p ±U q4 = 0. (2.6) 
x y 

Solving eqns. ( |2.5[) and ( |2.6| ) simultaneously we get the two mass ratios as 

x _ Ug?,Ub-iUp2Uq2 ~ Ug2Ub2UpiUqi ^ 7) 

U a lUblU p3 U q 3 — UazUbzUpiUqi + z(JJ a4 U b^P p zU q % — U a zUb1,U p4 U q4 ) ' 

y _ Ua3Ub3Up2Uq2 + ^a2^fc2^p3^3 & 2iS 13 ^ ^ 

UalUblU p 2U q 2 ~ U a2 Ub2U p iU q i + Z,{U a4 U b^U p2 U q2 ~ U a ?,Ub?,U p 4U q /l) 

The modulus of these quantities gives the magnitudes x m , y m while the argument deter- 
mines the Majorana phases a and (5. 

%m — |^| 5 Um — \y\ (2 - ^) 

a = arg(x), (3 = arg(y). (2.10) 

Thus, the number of the free parameters is five, the lowest mass mi (NH) or 7713 (IH), three 
Dirac and one Majorana type CP phases. We can check for the two mass spectra in terms 
of the magnitude of the mass ratios x m , y m and 

• SNH which corresponds to x m < 1 , y m < 1 and z m > 1 

• SIH which implies x m < 1 , y m > 1 and z m > 1 
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Thus, it is y m which determines if the hierarchy among the three light neutrinos is 
normal or inverted. Note that if the three light neutrinos are quasi-degenerate then we will 



have x m ph y m m 1. Unlike the three generation case discussed in [[D], 18] the lowest mass 
can not be determined in the four neutrino analysis in terms of x m and y m since these 
ratios also depend on m\ through z. Thus, we keep the lowest mass as a free parameter. 
To find out the allowed two zero textures we adopt the following procedure. 
We vary the lowest mass randomly from to 0.5 eV. All the five mixing angles in Table 2 
(apart from #34) and the three mass squared differences are distributed normally about the 
best-fit values with the la errors as given in Table 2. The three Dirac and one Majorana 
type CP phase as well as the remaining mixing angle #34 are randomly generated. Then, we 
use the above conditions to find out which mass spectrum is consistent with the particular 
texture zero structure under consideration. We also calculate the three mass squared 
difference ratios 

R v = 
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The 3cj ranges of these three ratios calculated from the experimental data are 

R u = (0.02 - 0.04), 
R ul = (1.98 x 10" 3 - 3.3 x 10" 3 ), 

R v2 = (0.63 x 10~ 4 - 1.023 x 10~ 4 ). (2.12) 

The allowed textures are selected by checking that they give the ratios within the above 
range. 



3. Results and Discussions 



In this section we present the results of our analysis. First we briefly discuss the results 
that we obtain for the two zero textures of the 3x3 mass matrices. Next we present the 
results that we obtain for the 3+1 scenario i.e 4 x 4 mass matrices. 

3.1 Results for 3 neutrino mass matrix 

For the 3 neutrino case, the lowest mass and the two Majorana phases can be determined 
from the mass ratios. Hence, the only unknown parameter is the Dirac type CP phase (^13) 
which is generated randomly. All the other oscillation parameters are distributed normally, 
peaked at the best-fit and taking their one sigma error as width. We find that all 7 textures 
which were allowed previously remain so. However, the textures belonging to A class allow 
NH whereas for the B class, B\ and B% admit NH and B2 and allow IH solutions. 
Class C gets allowed only for IH. The D, E and F classes remain disallowed. In the 2nd 



-7- 



y„oi 



0.001 0.01 ' 

Lowest Mass(eV) 




a( $oooi 0.001 o.oi o.i 
Lowest Mass(eV) 



.0001 0.001 0.0! 

Lowest Mass(eV) 



y„,(Ur 



0< $0001 0.001 0.01 0.1 

Lowest Mass(eV) 




y„o.i r 



0.001 0.01 

Lowest Mass(eV) 



y m o. 



0.001 0.01 

Lowest Mass(eV) 



y m o.i 



"'fi'oiKll 0.00! 0.01 0.1 

Lowest Mass(eV) 



01 0.001 0.01 0.1 

Lowest Mass(eV) 



y m o. 



0.001 0.01 0.1 

Lowest Mass(eV) 



Figure 2: The values of y m as a function of the lowest mass (mi or 7713) for the 3+1 case 
when the known oscillation parameters are varied in a Gaussian peaked at their 

respective best-fit values. 



column of Table || we summarize the results that we obtain for the two zero neutrino mass 
matrices with three active neutrinos using normal distribution of the oscillation parameters. 
The results obtained in this case are somewhat different from that obtained using random 
distribution of oscillation parameters. The reason for the difference stems from the different 
range of values of the atmospheric mixing angle #23 used by these methods. If we assume 
a Gaussian distribution for sin 2 623 around its best-fit then there is very less probability of 
getting the 3<r range in the higher octant as these values lie near the tail of the Gaussian 
distribution. This disallow B2 and B4 for NH and B\ and B3 for IH []2C|J. Similarly QD 
solutions for B class requires 623 ~ vr/4 JT6| ] and for a normal distribution of O23 with the 
peak at present best-fit the 3<r range extends upto ~ 44° and there is very little probability 
of getting values close to 7r/4. Similarly for the C class NH and QD solutions are allowed 
only for 823 values close to 7r/4 and hence is not admissible when Gaussian distribution of 
oscillation parameters about the best-fit value is assumed. 
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Table 3: The allowed mass spectra in 3 and 3+1 scenarios. The second and third columns give 
the results for the two zero neutrino mass matrices for three active neutrinos by assuming random 
and normal distribution of oscillation parameters respectively. The last column gives the allowed 
spectrum in the for the 3+1 case assuming normal distribution of parameters. For random distri- 
bution of oscillation parameters similar mass spectra get allowed although the parameter space is 
reduced in size. See text for details. 



3.2 Results for 3+1 scenario 

Adding one sterile neutrino, there exist in total forty five texture structures of the neutrino 
mass matrix which can have two zeros. 

(i) Among these the 9 cases with m ss = are disallowed as the mass matrix element 
contains the term m 4 U^ 4 which is large from the current data and suppresses the 



m 



other terms. Hence, m ss cannot vanish. 

(ii) There are 21 cases where one has at least one zero involving the mass matrix element 
of the sterile part i.e m ks = where k = e, fi, r. This element is of the form, 

m ks = miUkiUsi + m 2 U k2 U s2 e~ ia (3.1) 
+ m 3 U k3 U s3 e l ^~^ +m 4 U k4 U s4 e l(25 ^\ 

The last term in this expression contains the product m 4 U s4 which is quite large as 
compared to first three terms and thus, there can be no cancellations. Thus, neutrino 
mass matrices with one of the zeros in fourth row or column are not viable. 

(iii) The remaining cases are the 15 two zero cases for which none of the sterile components 
are zero. Thus, these also belong to the two zero textures of the three generation 
mass matrix. A general element in this category can be expressed as, 

m k i = miUkiUn + m 2 U k2 Ui 2 e~ ia (3.2) 
+ mzUksUtteWu-ft + m 4 U k4 U l4 e l ^~''\ 

here, k,l = e, /i, r. We find all these 15 textures, presented in Table [l] get allowed with 
the inclusion of the sterile neutrino. This can be attributed to additional cancellations 
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that the last term in eq. 3.2 induces. Table 3 displays the nature of the mass spectra 
that are admissible in the allowed textures. 

In Fig. 2 we present the values of y m vs the lowest mass for textures A±, B±, B%, C, 
D±, E±, E$, F\ and F%. This figure shows that for textures belonging to the A and E 
classes y m remains < 1. Thus, these classes admit only NH solutions. The textures 
belonging to the D class allow NH and IH while the B,C,and F classes allow NH, IH 
and QD mass spectra. 

The textures A\ — A2, B\ — B2, -B3 — -B4, D\ — D%, E\ — E 2 and F2 — F3 are related 
by P^ T symmetry where for the four neutrino framework P^ T can be expressed as, 

(l o\ 
10 
10 
\0 1 J 

in such a way that 

A 2 = t%A x P^. 

Note that for 3 generation case the angle #23 in the partner textures linked by fi — r 
symmetry was related as #23 = (§ — #23)- However, for the 3+1 case no such simple 
relations are obtained for the mixing angle #23- The angles 024 and #34 in the two 
textures related by fi — r symmetry are also different. For this case, the mixing angles 
for two textures linked by P^ T symmetry are related as 



p 

1 flT 



'12 



#12) #13 — #13> #14 — #14, (3.3) 



sin #24 = sin #34 cos #24, (3-4) 

. - COS #23 COS 6*34 - SU1 2 3 SU1 34 SU1 6> 2 4 fo c s 

sin 023 = , . (3-5) 

Jl - COS #24 sm #34 

• a sin6»24 ,„ c s 

sm^4 = 1 (3.6) 

Jl - COS ^24 sm #34 

The texture zero conditions together with the constraints imposed by the experimental 
data allow us to obtain correlations between various parameters specially the mixing angles 
of the 4 th neutrino with the other three for the A and E classes. For the B, C, D and F 
classes one gets constraints on the effective mass governing 0v/3/3. 

In order to gain some analytic insight into the results it is important to understand 
the mass scales involved in the problem. The solar mass scale is J Am^i ~ 0.009 eV 



whereas the atmospheric mass scale is yAm^ ~ 0.05 eV. Normal hierarchy among the 
active neutrinos implies m\ « m,2 « m% corresponding to m\ ^ 0.009 eV. It is also 
possible that mi ~ m 2 « m% implying m\ ~ 0.009 eV - 0.1 eV. We call this partial 
normal hierarchy. IH corresponds to m^ « mi ~ m 2 - If on the other hand mi > 0.1 eV 
then mi « 7712 ~ m% which corresponds to quasi-degenerate neutrinos. 
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Figure 3: Correlation plots for A and E class. 



• A and E Class 

For these classes we find y m to be mainly in the range > 0.0001 eV extending up to 
~ 0.1 eV. Thus, these classes allow normal hierarchy (full or partial) among the 3 
active neutrinos. These classes are characterized by the condition m ee = 0. m ee for 
the four neutrino framework can be expressed as, 

2 2 2 2 2 -ice 2 , 2 -iS 2 -i-y 2 /o rj\ 

mee = c 12 c 13 c M mi + c 13 c 14 e m 2 s 12 + c 14 e M m 3 s 13 + e 'm 4 s 14 . (3.7) 

For smaller values of m\ and NH the dominant contributi on to t he magnitude of the 
above term is expected to come from the last term s^Am^ ~ 0.022. Therefore, 
very small values of m\ is less likely to give m ee = for normal hierarchy. However, 
we get some allowed points in the small mi regime which implies smaller values of 
sf 4 . m ee can be approximated in the small mi limit as, 

m ee « e~ ta m 2 s\2 + e~ ll3 m 3 sl 3 + e'^m^sl^. (3.8) 

The maximum magnitude of the first two terms is ~ 0.003. Then using typical values 
of m.4 (~ 0.9 eV) from the 3a range, we obtain sf 4 ~ (0.003 — 0.004) in the small mi 
limit. This is true for all the textures in the A and E class. For the ^4i class we also 
simultaneously need m eA1 = 0. In the small mi limit approximate expression for m efl 
is 

m efl w ^""^-^14524^4 + e^-^s 13S23 m 3 + e- ia c 12 c 23 s 12 m 2 , (3.9) 

and the first term i.e m4Si4S24 ~ (0.05 — 0.06)s24- While the other terms are of the 
order (0.006 - 0.007) which implies s| 4 ~ (0.01 - 0.02). This is reflected in the first 
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and second panels of Fig. 3 where the correlation of sf 4 and S24 with mi is depicted. 
As mi increases the contribution from the first three terms in m ee increases and sf 4 
becomes larger for cancellation to occur. For m eA1 , this increase in si 4 helps to achieve 
cancellation for higher values of mi and therefore stays almost the same. Similar 
argument also apply to the Ei class which has m^ = 0. 

For A 2 class, in addition we have m eT = 0. In the limit of small mi, m er can be 
approximated as, 

m eT « e^-^-^suSum^ + e i ^-«s> 3 - rn 2 e~ ia si 2 ci 2 s 23 . (3.10) 

As discussed earlier m ee = implies small sf 4 ~ (0.002 — 0.005) in the limit of small 
mi. Thus, the contribution from the m 4 term is ~ (0.04 — 0.06)s3 4 . The typical 
contribution from the last two terms is ~ 0.008. This implies s| 4 to be in the range 
(0.02 - 0.04) for smaller values of mi. This is reflected in the third panels of Fig. 3 
where we have plotted the correlation of s| 4 with mi. Since with increasing mi, s i4 
increases to make m ee = 0, s| 4 does not increase further. Similar bounds on Sg 4 are 
also obtained for E 2 class. 

As one approaches the QD regime then the terms containing the active neutrino 
masses starts contributing more. So for higher values of mi complete cancellation 
leading to m ee = even at the highest value of sf 4 is not possible. This feature 
restricts mi to be < 0.1 eV in A class. 

Textures belonging to A and E class contains m ee = which is not possible for IH 
in the 3 generation case since the solar mixing angle is not maximal. In the 3+1 
scenario a cancellation leading to m ee = is possible for IH but it requires values of 
sf 4 in the higher side. It also contains a strong correlation in the Majorana phases. 
But in the other mass elements these constraints are not satisfied simultaneously and 
as a result the textures that contain m ee = do not admit inverted hierarchical mass 
spectrum. 

Since m ee = 0, the effective mass (m e // = |m ee |) governing the neutrinoless double 
beta decay (0v(3(3) is vanishing for these classes. 

• B, C Classes 

In the 3+1 scenario, B and C classes allow all three mass spectra - NH, IH and QD 
assuming the known oscillation parameters to be normally distributed (cf Table 3). 
In this case since #23 in the textures related by [i — r symmetry is not correlated in 
a simple way, the value of this angle not being in the higher octant does not play a 
significant role as in the 3 generation case. Among these only Bi allow few points 
for smaller values of mi for NH. For the IH solution, larger number of points are 
obtained corresponding to the lowest mass > 0.01 eV as is seen from Fig. 2. In these 
textures, for higher values of the lowest mass the active neutrino contribution to the 
matrix elements are larger and it is easier to obtain cancellations. Hence, textures 
belonging to these classes show a preference for QD solutions. For these textures the 
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Figure 4: The effective mass governing Qvj3j3 as a function of the lowest mass. The red 
(dark) points correspond to NH while the cyan (light) points correspond to IH. 



effective mass governing 0v(3f3 is non-zero. In the first row of Fig. 4 we present the 
effective mass as a function of the lowest mass for the textures B\, B% and C for both 
NH and IH. These two merge at higher values of the lowest mass corresponding to 
the QD solution. The effective mass in these textures is > 0.002 eV for NH and > 
0.02 eV for IH. If no signal is seen in future Qv(3f3 experiments then large part of the 
parameter space belonging to these textures can be disfavoured. 

• D,F Classes 

These two textures are disallowed in the 3 generation case. However for the 3+1 
scenario they get allowed. NH is admissible in all the textures belonging to these 
classes. The reason for this is the following. 

In the three active neutrino s cenario , the neutrino mass matrix in a fx — r block has 
the elements of the order of J Am^ « O.OleV for normal hierarchy. Thus, in general 
these elements are quite large and cannot vanish [p^ j. However, in the 3+1 case 
when there is one additional sterile neutrino, the neutrino mass matrix elements get 
contribution from the sterile part of the form m^UkiUn where k = e, fx, t. This term 
is almost of the same order of magnitude and thus can cancel the active part, resulting 
into the possibility of vanishing elements in the pt — T block. Thus, the zero textures 
which were disallowed for NH are now allowed by the inclusion of sterile neutrino 
(3+1 case). In the case of IH, m^ T element for three active neutrinos is always of 
the order of J Am^ ~ O.Oley and thus the textures D\, D2, F2, F3 which requires 
m iiT = were not allowed. However, for the 3+1 scenario the extra term coming due 
to the fourth state helps in additional cancellations and IH gets allowed in these (cf. 
Table 3). For F\ class, IH for three active neutrino is disfavoured because of phase 
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correlations. However, with the additional sterile neutrino this can be evaded making 
it allowed. In the bottom row of Fig. 4 we present the effective mass governing 
for the textures D\, F± and F2 as a function of the lowest mass. The texture D\ 
allows lower values of m\ for NH while for IH the lowest mass is largely k; 0.01 eV. 
QD solution is not allowed in D class . For F class more points are obtained in the 
QD regime. Future experiments on would be able to probe these regions of 

parameter space. 

The results presented above are obtained by varying the known oscillation parameters 
given in Table 2 as distributed normally around their best-fit and with a width given by 
the la range of the parameters. There is a finite probability of getting the points in the 
3a range of this Gaussian distribution although more points are selected near the best- 
fit values. However, note that for some of the parameters the 3a range obtained in this 
procedure is different from that presented in Table 2. Thus, our results may change if we 
vary the parameters randomly in their 3a range as we have seen in the 3 generation case. 
In Fig. 5 we show the allowed values of y m as a function of the lowest mass for the case 
where all the parameters are varied randomly in their 3a range. We find that lower values 
of the smallest mass get disfavoured by this method. The main reason for this is that if 
we use the Gaussian method then the allowed 3a range of the mixing angle sf 4 is from 
(0.002 - 0.048) while that of s| 4 is from (0.001 - 0.06). Thus, smaller values of sf 4 and 
S24 are possible which helps in achieving cancellation conditions for smaller values of mi 
or 771,3. But if the parameters are varied randomly in the 3a range presented in Table 2 
then such smaller values of the angles are not allowed and consequently no allowed points 
are obtained for smaller values of masses. In particular, we obtain mi(NH) or 771,3 (IH) > 
0.01 eV in all the textures However, main conclusions presented in Table [3| regarding the 
nature of the allowed mass spectrum for the 3+1 scenario remain unchanged though the 
allowed parameter space gets reduced. Specially for the A, E and C classes very few points 
get allowed. Fully hierarchical neutrinos (mi < 7712 < 7773) are not possible in any of the 
textures. Textures belonging to the B and F classes give more points in the QD regime. D 
class allows partial NH or IH. 

4. Conclusions 

Recent experimental observations make a case for enlarging the scope of three flavour 
oscillations to include one or more sterile neutrinos with mass around ~ 1 eV. Although 
induction of more than one sterile neutrino may provide a better fit to the oscillation data 
the cosmological observations may be more consistent with the three active and one sterile 
picture with the sterile neutrino being heavier. With the addition of one sterile neutrino 
the parameter space describing neutrino masses and mixing at low energy increases to 
include four independent masses, six mixing angles and six phases. The low energy mass 
matrix in the flavour basis now consists of 10 independent elements as opposed to six 
elements for the three generation case. It is well known that for the three generation case 
the neutrino mass matrix in flavour basis can have at the most two zeros. In this work, we 
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Figure 5: The values of y m as a function of the lowest mass when parameters are varied 

randomly. 



have considered the two zero mass matrices in the framework of three active and one sterile 
neutrino. We find many distinctive features in this case as compared to the three neutrino 
scenario. For the 3+1 case there can be 45 possible two zero textures as opposed to 15 
for the 3 generation case. Among these 45 possible two-zero textures only 15 survive the 
constraints from global oscillation data. Interestingly these 15 cases are the 15 two-zero 
textures that are possible for three active neutrino mass matrices. While for the three 
active neutrino case only 7 of these were allowed addition of one sterile neutrino make all 
15 cases allowed as the sterile contribution can be instrumental for additional cancellations 
leading to zeros. All the allowed textures admit NH. The classes B, C, F also allow IH and 
QD solutions in addition. The results are summarized in Table ||. 

If we vary the mass and mixing parameters normally peaked at the best-fit value and 
la error as the width then we find solutions for smaller values of mi(NH) and (IH). 
In this case for the textures with m ee = i.e A class and E class we obtain correlations 
between the mixing angles sin 2 #14, sin 2 #24 5 sin 2 #34 and the lowest mass scale mi. For these 
textures the effective mass responsible for neutrinoless double beta decay is zero. For the 
other allowed textures we present the effective mass measured in neutrinoless double beta 



decay as a function of the smallest mass scale. If however, the known oscillation parameters 
are varied randomly in their allowed 3a range then although the main conclusions deduced 
above regarding the allowed mass spectra in various textures remain the same the allowed 
parameter space reduces in size. In particular, we obtain a bound on the smallest mass as 
m sma iiest > 0.01 eV and completely hierarchical neutrinos are no longer allowed. 

In this work we have concentrated on the two zero textures. However it is possible 
that the neutrino mass matrix in the 3+1 picture may allow the presence of more than 
two zeros |^3|. These results can be useful in probing underlying flavour symmetries and 
also for obtaining textures of Yukawa matrices in presence of light sterile neutrinos should 
their existence be confirmed in future experiments. 
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